Photonic forces in the near field of statistically homogeneous fluctuating sources 
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Electromagnetic sources, as e.g. lasers, antennas, diffusers or thermal sources, produce a wavefield 
that interacts with objects to transfer them its momentum. We show that the photonic force exerted 
on a small particle in the near field of a planar statistically homogeneous fluctuating source uniquely 
depends and acts along the coordinate perpendicular to its surface. The gradient part of this force 
is contributed by only the evanescent components of the emitted field, its sign being opposite to 
that of the real part of the particle polarizability. The non- conservative force part is uniquely due 
to the propagating components, being repulsive and constant. Also, the source coherence length 
adds a degree of freedom since it largely affects these forces. The excitation of plasmons in the 
source surface drastically enhances the gradient force. Hence, partially coherent wavefields from 
fluctuating sources constitute new concepts for particle manipulation at the subwavelength scale 

PACS numbers: 42.50.Wk, 87.80. Cc, 42.25.Kb, 05.40.-a 



I. INTRODUCTION 

Photonic forces are increasingly studied due to their 
potential in many disciplines ranging from physics and 
chemistry to biology Of special importance is the 

manipulation of dipolar particles, understood as those 
for which the incident wave excites their first electric 
and/or magnetic Mie coefficients [HQ. Extensive stud- 
ies done on light from quasi-coherent sources show that it 
exerts mechanical action on these particles through both 
their conservative (gradient) and non-conservative com- 
ponents, allowing the design of optical tweezers which 
rely on the former component [2, 12, [|| & n d their recent 
extensions to the subwavelength , particularly nanomet- 
ric scale @, @|- On the other hand, the scattering, or 
radiation pressure, force component which until recently 
was believed to push objects d, 0, @-[H , has recently been 
designed to exert a pulling action towards the coherent 
source, as recently shown by exciting the induced mag- 
netic dipole or multipoles of the particle [T3, El , as well 
as by an appropriate design of the illuminating wavefield 
angular spectrum [l2| . 

We report here a new area of study for optical manip- 
ulation at the subwavelength scale, both theoretical and 
experimental, by partially coherent fields emanating from 
fluctuating sources [Lfl - fla ]. They completely change the 
nature of these forces and convey new behaviours to 
them. Interestingly, we find that planar sources, of such 
a general class as those that are statistically stationary 
and homogeneous, produce gradient forces that may be 
either attractive or repulsive. In turn, we demonstrate 
that these forces are dramatically enhanced as the co- 
herence length of the source decreases, as well as when 
surface plasmons (SPP) are excited on its surface. On 
the other hand, the non-conservative part of the force, 



composed of the radiation pressure plus spin density of 
angular momentum of the electric wavevector, is pushing 
and constant throughout the emission half-space. In this 
way, one can control the tractor or pushing effect of the 
resulting force on the particle according to the sign of the 
real part of its polarizability [12l. ll7H28| . 



II. FLUCTUATING OPTICAL FORCES 

Let us consider a fluctuating source emitting from the 
plane z = 0, (see Fig. 1). We shall assume that the ra- 
diated random field is described by an ensemble which 
is stationary, then we may work in the space-frequency 
domain 13[ so that its electric vector is expressed at fre- 
quency u) as an angular spectrum of plane waves propa- 
gating throughout the half-space z > \13L |29|: 



E(r,w) = / e(ks ± ,uj)e iks - r d 2 



S_L, 



(1) 



where k — uj/c, c being the speed of light in vacuum. The 
propagation vector k = ks is expressed as k = fc(sj_, s z ), 
so that sj_ = (s x ,s y ) are the transversal components of 
s and s z = — |s^| 2 , (|s^| 2 < 1), for homogeneous or 
propagating waves, and s z — i^\s\] 2 — 1, (|s^| 2 > 1) , 
for evanescent components . 

Let a dipolar particle with dynamic electric polarizabil- 
ity a e , be placed in the source vicinity. Being p = a e E 
the dipole moment induced in the particle by the E field, 
the ith Cartesian component (i = 1,2,3) of the mean 
force that the emitted wavefield exerts on it at frequency 
wis [li[13 
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FIG. 1. Illustrating the notation 



in the source plane z = is (3ljj £>y (pi, p 2) oS) 
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It is well-known [13| that Af k (ks±,ks' 



jk 

k 2 £jk{ks±,ks'±uj), where £jk (ks±, ks'±cu) is the four- 
dimensional inverse Fourier transform of (p\, p 2 ,a;). 
In addition, it was proven [l9| that for a homogeneous 
source the components of the electric angular correlation 
tensor are 5— correlated as 



A$ (ks ± ,ks' ± ,u) = k 4 S^ [k (s 



x£ 
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expressed as the sum of a conservative, or gradient force, 
pgrad p r0 p 0r ti nal ^ 5fta e and a non-conservative term 
i 7 !™ 2 proportional to 3a e . 5ft and 3 stand for real ad imag- 
inary parts, respectively. The symbol * denotes complex 
conjugate. The angular brackets mean ensemble average. 
Einstein's convention of omitting the sum symbol X)j=i 
on the repeated index j has been used. 

On introducing Eq. (p} into ([2]) one obtains 

Ff rad (r, oj) = -i~Ka e J f°° Tr^g (ks ±) ks' ± u) 

x (s* - s'i) e- lk ( s *- s '> r d 2 s ± d 2 s' ± , (3) 



<$( 2 ) representing the two-dimensional Dirac-delta func- 
tion. 

On introducing the above ^-function expression for 
(ks±, ks'±,w) into Eqs. ([3]) and (g]) one straight- 
forwardly obtains for the gradient force 

7 3 r 

= -i—^a e / H£<°) (fcs_L,o;) 

4 «/|s_lP>1 

x (s*-Si)e- 2k ^ 2 - u d 2 s ± , (6) 



ff c (r,w) = ^3a e 3 jiA ^ Tr^g (*»x, fca'^w) 

x ^e^ fc ( s *- s ')- r d 2 s± d 2 s' ± }, (4) 

(i,j,k = 1,2,3), Tr denotes the trace of the elec- 
tric angular correlation tensor .A*.^ (fesj_, ks'±,w) = 
(e*(fcsj_,o;)efc(A;s'j_,LLi)). Notice that since (E*E 3 ) is real 

and non-negative, Ff rad given by Eq. §3§ which equals 
jdta e di (E*Ej^ according to Eq. @, is a real quan- 
tity. Eqs. ([3]) and ((4]) reveal that whereas the gradient 
force depends on a weighted sum of the difference vectors 
s* — s' and, as we shall see, it has a negative sign if 5fta e 
is positive, thus pulling the particle towards the source, 
the non-conservative force associated to 3a e which is al- 
ways non-negative, only depends on the weighted sum of 
vectors s and pushes the particle forward. 

A. Statistically homogeneous sources. Gradient 
and non-conservative forces 

Let us address the wide variety of statistically homoge- 
neous sources, fl7t Il8l| . Then their electric cross-spectral 
density tensor" (ri,r 2 ,o;) = (E* (n) E } (r 2 )) 



The subindex in the integral of Eq. © means that the 
integration only extends to the non-radiative region be- 
cause the difference vector s* — s in Eq. ^ is clearly 
zero for propagating waves, (|s^| 2 < 1). Therefore 
the radiative components of the field emitted by statis- 
tically homogeneous sources do not contribute to the gra- 
dient force, which only depends on the evanescent com- 
ponents, ( | s_i_ | 2 > 1), for which s* — s = (0,0, s* — s z ) = 
(0, 0, — 2i-\/|s^| 2 — 1). Hence this force only exists in the 
near field, and depends on the distance z of the parti- 
cle to the source, having solely z- component normal to 
its surface. In addition, this force is attractive or repul- 
sive depending on the sign o/5fta e . Small particles with 
relative permittivity e > 1 have 5fta e > out of reso- 
nance and thus F^ rad (z,lj) will drag them towards the 
source. Conversely, near a resonance 5fta e may be nega- 
tive [3J, thus this force being repulsive. However, further 
study is required in this latter case, since then the par- 
ticle strongly scatterers the field emitted by the source, 
and therefore the analysis developed here should not be 
exact due to multiple scattering of the radiation between 
the source and the particle. Hence it is shown that the 
gradient force near a statistically homogeneous source is 
entirely of non-radiative nature and may work as a trac- 
tor force 

Analogously, from Eq. (jj) one also derives for the non- 
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conservative force F™ c a dependence on z only: 
2 
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F™£ and denote propagating and evanescent wave 

contributions, which correspond to the first and second 
integral terms of Eq. (0 , respectively. Notice that F"^ > 
is constant throughout z > . 

Let the source also be statistically isotropic [13J so that 



£|, 0) (Pi,P2,w) = £\"> {p,u>), where p = |pi - p 2 \- The 



,(0) 



spatial coherence function of the field in z = is [3lM3 

Tr£ jj (p, lj) and the spectral degree of spatial coherence 

p^ ( p ,lj) = Tr£^ } ( p ,lj) /SW(w), where the wavefield 

spectrum on the source is: S^°\oj) — TrS*?' (0, w). 

To illustrate these results, we shall consider a Gaussian 
spectral degree of coherence p^ (p, u>) — exp [— p 2 /2a 2 ] , 
so that taking Fourier inverse one obtains 



Tr£$ (ks x ,uj) = S^(uj)p^(ks ± ,oj) 



= S<°>(w) (a 2 /2n) 



exp 



-(fca|s ± |) 2 /2 



(8) 



where a is the correlation or coherence length of 
the source. On introducing this expression for 
Tr£-^ (fcsj_,u;) into the force equations ((SJ) and ©, we 
obtain that on writing in cylindrical coordinates: s x = 
sj_ cos 4>, s y — s±_ sin <fr, and due to the rotational sym- 
metry of the source, the transversal components of the 
non-conservative force are zero, viz. : F™ c y (z,uj) = 
since so are the corresponding integrals of Eq. ([7]) when 
one performs the azimuthal angle <fi integration. Also, 
since s z — i-^\sj_\ 2 — 1 for |s^| 2 > 1, the second inte- 
gral in Eq. ([7J is purely imaginary, which implies that 



z,ev 



0. 



is due to the non-radiative plane wave components of the 
emitted field, whereas to the latter only the radiative com- 
ponents contribute. As the distance from the particle to 
the source plane grows to values z > X, F^^f (z, ui) tends 
to zero due to its evanescent wave composition. Never- 
theless, as we shall see, the source coherence length a 
plays an important role on these contributions. 

The integration of Eqs. © and © using the Gaussian 
spectral degree of coherence, quoted before: p^ (p, ui) — 
exp [— p 2 /2a 2 ] , leads to an analytical expression for the 
gradient and for the non-conservative force. For the lat- 
ter, Eq. ([T7l) yields the proportion of radiation pres- 
sure and curl components for unpolarized emission. This 
calculation is straightforwardly done on making: s x = 
sj_ cos <j), Sy = sj_ sin <f>, and leads to 



F3 rad (z,w) = MaeS^^e-^ 2 * 2 ^ 



72^ 

V _3 



2(7 



)e^erfc(V2z/o-)}. 



(10) 



F r z lc {z,uj)=^a e S { -°\ui)/2[k 
1 



a V 2 



crfi(fc(T 



(11) 



where erfc(x) = 1— erf(x), erf (a;) being the error function: 
erf(a;) = 2/^/tt Jq e~ f dt, and erfi(x) is a positive real 
function defined as erfi(x) = evi(ix)/i. 



B. The curl force 

It is well-known [H, [3(| that the non-conservative part 
of the force F^ c is the sum of a scattering force, or radi- 
ation pressure 

F? c = (fc/2)3a e 5R(E x B*). ( 

= (l/2)9fae9f {(E^diEj) - (E*d 3 E t )} , (12) 

given by the averaged field Poynting vector, plus the curl 
of a electric spin density: 



Y 




Ti-£)° k > (ks ± ,uj) Sl d 2 s ± } . (9) 



Thus, while F™£ (z,to) > is constant throughout z > 0, 

as so is the spectrum S^(oj) propagating into z > 
[22| . the evanescent waves do not contribute to the non- 
conservative force F™ c (r,w). 

In summary there are therefore two force components 
acting on the particle: F^ r ^f (z, u) and F™ c h {z , u) , per- 
fectly distinguishable from each other since the former 



F nc,cun = (^2)3^3 ((E* • V) E)< 

= (l/2)9fa e 9f(^a,^>. 



(13) 



If the field emitted by the source is unpolarized: 
£^(p,Lu) — F^°\p,cj)Sjk, F^ oS> (p,Lu) being a scalar spa- 
tial correlation function whose two-dimensional Fourier 
transform will be denoted as F^°'(ks±, ui). Then 



Tr£$(ks ± ,u>) = 3FW(ks ± ,Lv), 



(14) 



and the radiation pressure contribution F™ c ' pr to the 
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non-conservative force is: 

pnc,pr 



s_l|<1 



Sj_|<1 



TvEjk (fcsj_,w) Si - £ji (ks±,uj) Sj 



(ks±) - (fasj.) 
k 3 %a e I F (0) (jfcsx) = F. 
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s»crs 



nc,pr 
z ' 



d 2 s 



(15) 



since the azimuthal angle integrations when s% is either 
s x or s y is zero. 

In a similar manner, the curl force contribution 

pnccurl tQ pnc ig 
1,3 

Fine, curl ^ rv 
i = 

fcV 
= y 9a e 



S I <1 



f ji (fcsj_, w) Sjd 2 s± 



f F^ (ks ± )s z d 2 s ± =F r z 

J\s ± \<l 



Namely, for unpolarizad radiation: 



pnc,pr t^pnc.curl JT nc 



(16) 



(17) 



III. EXCITATION OF SURFACE PLASMON 
POLARITONS. NUMERICAL RESULTS 

Without loss of generality, we shall also address surface 
plasmon polaritons (SPPs), excited on the source plane 
z = 0. Let this be gold for example, choosing for instance 
A = 459. 9nm, its permittivity is e = -2.546 + i3.37 [H. 
The SPP wave vector ks ± = fcs^ pp = ±k [s/ (e + 1)] 1/2 
corresponds to a pole of the Fresnel coefficient, (either 
on reflection or on transmission depending on the set-up 
configuration), R(ks±,uj) [29|, [35|. Then it is easy to ob- 
tain that the former equations ([6]) and (JT]) are valid on 
substituting £ ^ (ksx, u>) by £r£ (ks±, uj) \R (ks±, u>)\ 2 

M 

Figure 2 shows the attractive gradient optical force 
due to evanescent components for two random sources: 
one without and one with excited SPPs (cf. Fig. 2a 
and Fig. 2b, respectively). The normalized value 
Ff arf (z,w)/(fcS*(°) (cj)5Ra e /27r) is represented in arbi- 
trary units. As predicted by Eq. the gradient force 
drags the particle towards the source plane; (notice that 
since this normalization does not include 5ia e , it does 
not contain an eventual negative value of this quantity). 
In both figures we observe its exponential increase as the 
distance z of the particle to the source decreases. Nev- 
ertheless, this force is mainly governed by the coherence 
length a. For a — A/8 (red line), the magnitude of this 
force is maximum, but we observe that around a = A/2 
(blue line), and beyond, there is an important decrease, 
with values between 10~ 2 and 10 -3 in the magnitude 
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FIG. 2. Pulling gradient optical force due to evanescent com- 
ponents, a, Gradient force in arbitrary units (a.u.) versus 
distance to the source z/X for different values of the source 
coherence length a. b, The same force when SPPs are excited 
in the source. A significant decrease of the magnitude of this 
force is clearly seen as a grows about a — A/2. 



even at subwavelength distances z of this force, which is 
practically zero, (F~/ rad ~ 10~ 8 ), for a = A and z = 
(this latter curve is not shown). With this, we demon- 
strate that the decrease of the source coherence length 
gives rise to an increase of the gradient force and its ef- 
fect is larger than that of the distance z of the particle 
to the source plane. Eventually, a <5-correlatcd source, 
(then a — > 0), like e.g. a thermal source, will max- 
imize this force. In addition, we show with Fig. 2b 
that the excitation of SPPs in the source increases the 
strength of this near field force by approximately one or- 
der of magnitude. This is due to the then larger values 
of £p) (ks±,uj) \R (fcsj_, lu)\ 2 stemming from the pole of 
\R{ks ± ,io)\ 2 at ks s ± pp . 

Correspondingly, Figures 3a and 3b show the nor- 
malized total force F* ot (z,uj) = (2n/kS^ (ui)) ■ 
(Fs rad (z,uj)/^a e + F z nc (z,cj)/3?Q e ), in arbitrary units, 
without and with SPP excitation, respectively. At large 
distances (z > A), the total force is a constant of the 
distance z and repulsive according to the behaviour of 
the non-conservative component F z c , which dominates 
in this region of z, regardless of the value of a. In addi- 
tion, this non-conservative force is maximum for a = A/2, 
in contrast with the decrease of the gradient component 
as a increases. 

One might think that, due to its evanescent wave com- 
position, the magnitude of the gradient force at sub- 
wavelength distances would be higher than that of the 
non-conservative force, however this is not totally truth 
due the larger effect of the source coherence length on 
Fl rad rather than on F r z lc . In near-field F tot ~ F grad for 
°~ < A/4; but as a increases, F grad becomes negligible, 
being for a > A/4 F tot ~ F nc . These effects appear in 
Figs. 3a and 3b. Particulary, we see in Fig. 3b that 
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FIG. 3. Optical total force, a, Total force in arbitrary units 
(a.u.) versus distance z/X to the source in arbitrary units for 
different values of the coherence length a. b, The same as in 
(a) when SPPs are excited. In this second case we observe an 
increment Az at which the magnitude of the gradient force 
starts to exponentially increase 



IV. CONCLUSIONS 

We have reported a new area of study on photonic 
forces exerted on small particles by discussing near field 
effects due to fluctuating sources. This opens new per- 
spectives on subwavelength effects and manipulation in 
such general physics cases that range from light propa- 
gation through the turbulent atmosphere [36] , to speckle 
patterns from a large variety of statistical sources, also 
including scatterers, optical diffusers (37rl39| . as well as 
thermal or blackbody sources, which opens new possi- 
bilities at the subwavelength scale, particularly at the 
nanoscale. We have seen that in the large variety of sta- 
tistically stationary and homogeneous sources, only the 
evanescent components contribute to the gradient forces, 
while the non-conservative part that contains radiation 
pressure and curl forces is due solely to emitted propagat- 
ing components. Hence the subwavelength information is 
encoded in the gradient forces. Same numerical examples 
were given for statistically isotropic unpolarized emitted 
wavefields, showing the important effect that the source 
coherence length has on these forces, specially on the gra- 
dient component. Also the excitation of surface waves 
importantly enhances the magnitude of these forces. We 
expect that these findings stimulate experiments and ap- 
plications by this novel particle manipulation scenario. 
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